The Stability of the Low-dimensional Mixtures of Dilute Quantum Gases 
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We study the stability of the dilute Bose-Fermi and Bose-Bose mixtures with repulsive interactions 
in one and two dimensions in terms of the renormalization group. For the Bose-Fermi mixture, we 
show that the uniform mixture is stable against de-mixing in the dilute limit. For the Bose-Bose 
mixture, we give the stability conditions in the dilute limit. As a byproduct, we also calculate 
the critical temperature for the superfluid phase of the two-dimensional Bose-Fermi mixture in the 
extremely dilute limit. 
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I. INTRODUCTION 

Ever since the achievement of Bose-Einstein condensa- 
tion in the dilute gases of ultracold atoms, the study of 
degenerate quantum gases has become a subject attract- 
ing a lot of experimental and theoretical vigor—. Due 
to the experimental facilities which allow control of var- 
ious parameters characterizing the system, it provides 
an ideal laboratory for the studies of quantum many- 
body physics. Among them, one line of researches is to 
study the mixture of quantum gases. These include the 
mixture of bosonic and fcrmionic atoms (the Bose-Fermi 
mixture) 2 -, and the mixture of two different hyperfine 
states of the same bosonic atoms^ or two different kinds 
of bosonic atoms^ (the Bose-Bose mixture). The study 
of the degenerate mixtures of quantum gases at lower di- 
mensions is particularly interest and a lot of theoretical 
investigations are devoted to it£~— . This is because the 
competition between the inter-species and intra-species 
interactions and the strong quantum fluctuations in low 
dimensions may introduce many interesting phenomena 
such as phase transition, new quantum states, and quan- 
tum phase transitions. 

The present work is to study the stability of the low- 
dimensional Bose-Fermi (BF) and Bose-Bose (BB) mix- 
tures with repulsive interactions in the dilute limit. For 
binary mixtures with mutual repulsions, the system may 
experience a de- mixing transition by varying the strength 
of the interspecies interactions. The de-mixing transi- 
tion or the stability of the one-dimensional (ID) uniform 
mixture has been analyzed by using cither the mean-field 
theory^ or the bosonizationSii. The mean- field theory has 
the well-known problem that it can not treat properly the 
strong fluctuations in low dimensions. Though there is 
nothing wrong with the bosonization approach in one di- 
mension, a naive extrapolation of the results obtained 
by the bosonization to the strong-coupling regime may 
sometimes lead to incorrect conclusions. For example, 
the bosonization approach predicts that the BF mixture 
will be phase separated in the strong repulsion regime^. 
In the solvable limit, however, the Bethe-ansatz result 
indicates that the uniform BF mixture is stable against 
de-mixing irrespective of the strength of repulsions^. 
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FIG. 1: (color online) The schematic phase diagram of the ID 
(2D) BB mixture with repulsions. We have taken ni = n/2 = 
ri2- The dashed line is the mean-field value. In ID, the critical 
strength Uc at low density for mi 7^ ni2 is a linear function 
of n. In 2D, the critical strength wi 2 ' at extremely low den- 
sity for mi 7^ 1712 is a linear function of 1/ln [l/^na 2 ,)], where 



For mi 



m.2 



both U 



(i) 



and wi 2 ' take 



0. Moreover, and will 



m 2 / m l 
" — 2n(l + m 2 /m 1 ) 2 

a constant value u c as n 
approach the mean-field value u m f = ^/«nW22 at large den- 
sity. In between, we assume that they are smooth functions. 
In the SU(2) symmetrical limit, u c — u m f ■ 



Therefore, it is desirable to study this issue from a dif- 
ferent approach. For the ID single-component dilute 
Bose gas, it has been shown that the low-energy physics 
in the strong coupling regime can be well-captured by 
the renormalization-group (RG) method starting from a 
zero-density quantum critical point (QCP)^. Here we 
extend this idea to the low-dimensional mixtures with re- 
pulsions. Our results should provide a point of view com- 
plementary to the bosonization approach in the strong- 
coupling limit. 

In one dimension, there is an interacting zero-density 
fixed point which controls the low-energy physics of the 
ID uniform mixture in the strongly repulsive limit. Simi- 
lar to the single-component dilute Bose gas, this nontriv- 
ial fixed point can be captured by the e expansion, where 
e = 2— d. For the two-dimensional (2D) mixtures, all cou- 
plings are marginally irrelevant in the sense of RG. Thus, 
we expect that the stability conditions given by the mean- 
field theory are correct to the leading order. Our main 
results are as follows: (i) We find that the ID uniform 
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BF mixture in the dilute limit, i.e. the strong-coupling 
regime, is always stable against de-mixing for any finite 
strength of repulsive interactions. This is consistent with 
Bethe-ansatz result in the solvable limit^. And our anal- 
ysis extends this result to the region beyond the solvable 
limit, (ii) We derive the stability condition of the ID 
uniform BB mixture in the strong-coupling regime (Fig. 
[T]), which is different from the bosonization result, (iii) 
In the extremely dilute limit, we find that the interac- 
tions for the 2D mixtures are strongly renormalized so 
that the stability conditions obtained by the mean-field 
theory are changed considerably, (iv) We also calculate 
the critical temperature for the 2D BF mixture in the 
extremely dilute limit [Eq. (|26|) ] to investigate one as- 
pect of the effects arising from the inclusion of fermions. 
Although the stability of the mixtures with attractions is 
also an interesting issued, it is beyond the scope of the 
present work. This is because the fixed-point structure 
for attractions is different from the one for repulsions, 
i.e. the former exhibits a runaway RG flow, and thus our 
approach cannot give definite results. Finally, we notice 
that a similar RG analysis of the stability of the dilute 
BB mixture in low dimensions was also carried out in 
a recent paper—. The results arc different from ours in 
certain aspects. We will address the origin of this dis- 
crepancy at the end of the paper. 

The rest of the paper is devoted to the detailed anal- 
ysis that leads to the above conclusions. Sec. |TT] and 
Mil are about the BF and BB mixtures, respectively. We 
compare our results with the previous work in the last 
section. An appendix about the mean-field treatment of 
the BF mixture in d dimensions is provided for reference. 



II. THE BOSE-FERMI MIXTURE 

To study the mixture of bosons and spinless fermions 
in the dilute limit, we start with the action 



dr j d d x^(d T - ^— - fi b j <j) 
+ jy T Jd d x(?f\^+g bf \cbm 2 



assume that 



(1) 



where /3 = 1/T. (We will take h = 1 and k B = 1.) 
Here </> and tp are the fields describing the bosonic and 
fermionic atoms, respectively, m/, fif and m&, fib are 
the masses and chemical potentials for the bosons and 
fermions, respectively. The couplings gbf and gbb are 
related to the scattering lengths a&/ and abb through 
the relations gbf = 2uj±abf and gbb = 2oj±abb for the 
ID trap^ and gbf = y '4n(mb + rrif)uj±/nibTnfabf and 
gbb = \^omJ±Jrribabb for a 2D trap^ where ui± is the 
transverse confining frequency. We shall take gbb-, gbf > 
to avoid the possible formation of bound states. We also 



fib, \if < 



m a al 



(2) 



where a = b, f and ao is the average range of interac- 
tions between atoms, which can be regarded as the short- 
distance cutoff of the action S. We will see later that this 
condition amounts to the dilute limit. 

The action S [Eq. Q] has a QCP at fi = 0, T = 0, and 
9bb = = gbf (the Gaussian fixed point). In the dilute 
limit, this Gaussian fixed point may be an appropriate 
departure point to study the low-energy physics of this 
system. To proceed, we shall employ the RG method. 
The strategy is as the following. We first analyze the 
fixed-point structure of the theory. Then we integrate the 
RG equations to the scale where the correlation length 
becomes of 0(1) in the unit of the cutoff of the resulting 
renormalized Hamiltonain. At this scale the quantum 
fluctuations have been taken into account (within the e 
expansion) , and thus we may apply the mean- field theory 
to the renormalized Hamiltonian to study the stability 
of the mean-field solution corresponding to the uniform 
mixture. In terms of the solutions of RG equations, we 
may express the stability conditions in terms of the phys- 
ical quantities. 



A. Renormalization-group theory 



To proceed, we first make a change of variables: 



X = UqX 

-d/2 



t = ma n T 



a 



o 



-d/2- 



such that x, f , if), </> all become dimensionless, where m 
is a quantity carrying the dimension of mass. Then, the 
action S [Eq. ([1])] at T = can be written as 



S = jdTjd d i^[d f 

+ [df[d d xft( d f - 



2m/ 
V 2 



Tf U 



2rhb 



df d 



\bb 



2K d m b 



- n 

K d m b f 



where my = 2nibmf / (nib + m f), TO fc = rrib/m, fhf = 
nif/m, fhbf = m-bf/m, and 

r f = mfifal , r h = mfi b al , 
Afc6 = K d m b gbbao , Xbf = K d m b fgbfa% , (3) 

arc all dimensionless parameters. In the above, we define 
K d = 2/[(4 7 r)' i / 2 r(d/2)] (T(x) is the Gamma function) 
and e = 2 — d. We will see later that m always appear 
in the physical quantities in the guise of ma\. There- 
fore, for universal quantities which do not depend on a 
explicitly, the dependence on m is completely dropped 
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out. The determination of the short-distance energy scale 
1/(tooo) relies on the comparison with exact solutions or 
experimental results. 

To perform the RG transformations, we decompose the 
fields ip and </> into tj) = ^< + ^> an d 4> — 4>< + 4>> \ where 
ip > ((f> > ) consists of the Fourier components of ip (<p) with 
Ae~ l < \k\ < A. and ip< (</><) consists of the Fourier 
components of ip (cj)) with \k\ < Ae~ l . Here A is the 
(dimcnsionless) UV cutoff for the momentum and I > 
is the scaling parameter. By integrating out ip > and </>> 
to the one-loop order and then rescaling the momentum 
and frequency by k — > e~ l k and tu — > e~ 2l u>, we obtain 
an effective action for the slow modes. By ignoring the 
irrelevant operators in the sense of RG, this effective ac- 
tion has an identical form to the original action S but 
with different values of the parameters r a (l), Xbb(l), and 
A{,/(7), which can be regarded as the corresponding pa- 
rameters at the momentum scale Ae~ l . By comparing 
the effective action of the slow modes with the original 
action and considering an infinitesimal value of Z, we get 
the one-loop RG equations 



(4) 



and 



cZA 



hi, 



dl 
d\bf 



— eXbb — 



\2 

A bb 



1 - 2m b r b 



X 2 
A bf 



1 - fribf{rb +rj) 



(5) 



where the initial values r Q (0), At>b(0), and Ab/(0) are 
given by Eq. ([3]). Since we are only interested in the 
dilute limit so that r a (l) <C 1 before scaling stops, Eq. 
([5]) can be approximated as^ 



(6) 



where A represents either Abf, or Xbf . 

If d ^ 2, Eqs. (H]) and ^ will have two fixed points: (i) 
the Gaussian fixed point (r Q , Xbb, Xbf) = (0, 0, 0) and (ii) 
the interacting zero-density fixed point (r a , Xbb, Xbf) = 
(0, A*, A*) with A„ = e. On the other hand, for d = 2, 
the two fixed points merge into one - the Gaussian fixed 
point. Since gbb,9bf > 0, in the dilute limit, the zero- 
density fixed point will control the low energy physics 
in one dimension, whereas in two dimensions it is the 
Gaussian fixed point which will govern the low energy 
physics. Since the low density implies strong couplings 
in one dimension, this nontrivial zero-density fixed point 
will describe the strong-coupling physics. 

The solution for Eq. ([6]) in d < 2 is of the form 



A»(l) 



A* 



1 + CbbC 



with 

Cbb 



A* 



K d m b g bb al 



1 , c bf 



; - A * _ d , (7) 
1 + Cb/e 61 



A* 



while in d = 2 it takes the form 
K 2 m b gbb 



X bb (l) 



1 + K 2 m b gbb 



x m _ K 2 m bf g bf 

- , Xbf(i) - j ■ (8) 

/ 1 + K 2 m b fgbfi 



Scaling stops at I = Z* where Eq. ([6]) can not be applied 
anymore. This occurs when rb(l*)+rf(h) = 1,— yielding 



(9) 



From Eq. we see that the condition ^ implies that 
e l * > 1. 

The uniform mixture corresponds to the state with 
nonvanishing values of p b = (|<K^*)| 2 ) an d Pf = 
(\4>{l*) | 2 )- Within the framework of the RG, the phase 
boundary of the uniform mixture is determined as fol- 
lows: First of all, we calculate the effective potential of 
the order parameter pb at given \i a : 



-{d+2)l, _ 

f = — zpr/O. 



ma 



(10) 



where f(l*) is the effective potential of the renormalized 
Hamiltonian at / = Z*. Next, as a function of pb, the equi- 
librium state is determined by the absolute minimum of / 
at fixed [i a . Since the prefactor in Eq. (fT(j|) is always pos- 
itive and does not depend on pb, the absolute minimum of 
/ corresponds to the absolute minimum of /(Z*). There- 
fore, to search for the absolute minimum of /, it suffices 
to consider /(Z*). In particular, the stability of the uni- 
form mixture can be determined from /(/*). Because the 
coupling constants at Z = Z* arc small (in the sense of the 
e expansion for d < 2), we may calculate /(Z») in terms 
of the perturbation theory in A(Z»). The leading order 
result will given by the mean-field theory. In the follow- 
ing, we will employ this RG-improved mean-field theory 
to study the stability of the uniform mixture. To avoid 
the confusion, we emphasize again that / is a functional 
of the order parameter pb at given ^, a , and its absolute 
minimum determines the thermodynamical relation be- 
tween pb and fi a and the associated stability conditions. 
The free energy density fi is obtained by inserting the 
resulting value of pb into /, i.e. £1 = f(n a , Pb(Ha))- 



B. The uniform mixture in one dimension 

We first study the stability of the ID uniform mixture. 
From Eqs. (|A2j) and (|A3|) . to the leading order, pb and 
Pf are determined by the mean-field equations 

^bf * Xl b 

js ~ PS = r b- > 



and 



K d(„~ * {™b + mf)Xl f 
p f = -{ 2m ^ K^b— Rb 



d/2 



4 



where r* = r Q (7»), A£ b = X b b(h), and A 
In the dilute limit, we may simply set A 



bf 



hi, 



A* 



Xbf(l*). 



X 



bf 



Hence, by solving the two equations within the e expan- 
sion, we get 



Pb 
Pf 



K d rh b 

X* 
2K d rhf 



I - ~~r(™* + rh f) r } + °( e ) ) 

r} - —f(™>b + rhfVl + 0(e) . (11) 



We still have to relate fj, a (or r b and rj) to the bo- 
son density rib and the fermion density n/. This can be 
achieved through the scaling equations 



n b 



-da. 



_ -dh 



e~ dl 'a^ d p b , 
-°- dU % d p f 



(12) 



Inserting Eqs. ([9]) and (|11|) into Eq. (fT2|) and noticing 
(if), we find that 



that r*= fi a /(fib 
X 



j. 



and 



2K d m f 



d(d + 2] 



AK d fhf 

2£fl6+(l + 0»/ 



(13) 



(l + 3£)n 6 + [d(d + 2)/(2A*)]n/ 
(1 + Qnb + [d(rf + 2)/(4A,)]n/ 
(1 + 3£)n 6 + [d(d + 2)/(2A,)]n/ 



(14) 



where £ = mf/mb. By comparing Eq. (|9|) with Eq. 
(|13[). we see that the condition <j2j) indeed amounts to 
the diluteness condition n Q aQ -C 1. 

Now we are able to study the stability condition for 
the ID uniform mixture. From Eq. (|A5[) . the stability 
condition can be written as 



e/d _ 

pf > 



Kd 
d 



^ d m b m f (X* bf f 



l bf^bb 



Inserting Eqs. 0, (jT2]) . and ((T3j) into Eq. 
stability condition can be written as 



1 > X. t G(n b /rif) 



c bb e 



(l + c &/ e- d -) 2 



(15) 
the 

(16) 



where 



G(x) 



(i + O 2 

4£ 



d + 2 



We may use Eq. (|16jl to determine the critical strength 

9hf( e ) °f 9b f at given g bb , nb, and rif by assuming that 
eCl. Then, the critical strength of g b f in e? = 1 can be 
estimated by setting e = 1. To proceed, we notice that 
e -eZ* mus ^ treated as an exponentially small quantity 
which can not be expanded as a power scries in e. This is 
because Eq. (|16|) arc obtained by the double expansion 
in e~ d * and £j21 Moreover, G(n b /nf),c bb = 0(1) <C 1/e 



for given c/^f,, ^b, and n/. Keeping these in mind, we find 
that 



K d m b fg b fa e > 



-eh 



^eG(n b /n f ) 



(17) 



Since the value of the R.H.S. in the inequality (fTT|) is neg- 
ative, it is always satisfied as long as gbf > 0. This means 
that the uniform mixture in 2 — e dimensions is always 
stable for any finite strength of g b f > 0. By extrapolat- 
ing this result to e = 1, our analysis suggests that the 
ID uniform mixture may be stable against de-mixing for 
any finite strength of repulsions. 

The action S [Eq. (JTJ] can be exactly solved by 
the Bethe ansatz under the conditions m b = rrif and 
9bb = 9bj s - In the solvable limit, the uniform BF mix- 
ture is shown to be stable against de-mixing. This result 
is correctly reproduced in our approach, and is extended 
to the region beyond the solvable limit. 



C. The uniform mixture in two dimensions 

Now we turn to the 2D case. Here we will study two 
problems: (i) the stability condition for the uniform mix- 
ture and (ii) the effect of fcrmions on the critical temper- 
ature T c of the superfluid (SF) phase. 

We first examine the stability condition. Following 
from Eq. (|T5|) . the stability condition in d = 2 becomes 



(K f ) 2 < 



2m bf 
m b + ni 



■X 



bb 



(18) 



The diluteness condition only requires that e l * 3> 1. On 
account of the logarithmic RG flow for the coupling con- 
stants in d = 2, we have two situations: 

When e'* ^> 1 but Z» = 0(1), we may set the cou- 
pling constants by their bare values in Eq. (|18|) . i.e. 
X bb = K2m b g b b and A^ = Kim b fg b f. Thus, the stability 
condition in this case is 



2 . 27r .9bfc 
9bf < 



(19) 



This is just the mean-field result. 

On the other hand, when Z* 3> 1, which is referred to 
as the extremely dilute limit, we have to insert Eq. ([8"j) 
into Eq. (fig]), yielding 



L > 



(m b + mj) s 
4m b m f 



1 - 



8tt 



(m b + m f )g b f 



(20) 



We still have to determine Z* as a function of n a . For 
Z* ^> 1, we may simply set A£ b = 1/Z* = X b j. The rest of 
the procedure is similar to the ID case, and the result is 



-ml 

2 V 27m / a o 
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to the leading order in the 1//* expansion. Consequently, 
the extremely dilute limit is given by the condition 



lnl 



1 



( 

\2irnfaQ 



> 1 . 



(21) 



Since the right hand side of Eq. (|20|) is of 0(1), it is 
automatically satisfied as long as Eq. (|2"Tj) is obeyed. 
Therefore, we conclude that in the extremely dilute limit 
the uniform mixture is stable against demixing for any 
finite strength of repulsions. In other words, the quantum 
fluctuations stabilize the uniform mixture with repulsive 
interactions in the extremely dilute limit. 

In the parameter space where the uniform mixture is 
stable, the system will becomes a supcrfluid (SF) phase 
through the KT-typc transition as lowering the tempera- 
ture T. Here we would like to calculate the critical tem- 
perature T c for the SF phase in terms of the RG. Now 
there are three relevant parameters r&(Z), rf(l), and i; in 
the sense of RG, where the dimensionless temperature ti 
is defined as t\ = toe 21 with to = mTa\. To proceed, we 
will assume that \n a \ T c <C l/(m Q a,g). In this param- 
eter regime, it is U which reaches 0(1) first under the RG 
transformations. Hence, we run the RG equations to the 
scale I = I such that tj = 1. Accordingly, I is given by 



Here, we keep only the leading term in the small f a ex- 
pansion. Inserting this result into Eq. (|23jl . we obtain 



h w -2Aln(2A) 



(24) 



We have to express f& as a function of n a and T. The 
procedure is similar to the ID case, and the result is 



A 



A 



n 



K 2 mbT KirribfT 
Inserting this expression and using Eq 



we find 



m b T \ 



2* 



-71/ 





1 




) =ln {l ln 


} 




mTa 2 , 





From Eq. (J25J), T c is given by 

7r[TH+(l + 0/(20n/] 



m b ln{± 



hi 



(25) 



(26) 



By setting n/ = 0, Eq. (f2"6")l reduces to the critical tem- 
perature for a single-component dilute Bose gas^ Equa- 
tion (|26|) suggests that the inclusion of fcrmions increases 
the critical temperature. 



I 



In 



mTa 2 ) 



(22) 



Before I reaches I, ti <C 1, and thus we may approximate 
the RG equations for r a (l), Xbb(l), and Xbf(l) as the ones 
at T = 0. This may introduce multiplicative errors of 
order unity coming from the imprecise treatment of the 
regime ti ~ 1, but makes the analysis simpler. At I = I, 
the action S becomes 



df / d 2 x if}^ 



df 



2m 4 



1> 



df d 



> 2 x 



df I d 2 £ 



V 2 

Of - — n 

2m b 



X 1114 1 A -|0| 2 |^| 2 



2K 2 m b K 2 m b f 



where f a = r a (l) = [L a jT and A = I /I. Here we assume 
that I ^> 1 so that X b b(l) ~ A and \f(l) ~ A. Since 
A< 1, a perturbativc expansion in A is reliable. 

The field <f> may be considered the order parameter of 
the SF phase. Hence, the critical temperature can be 
determined by the equation 



h + E = 



(23) 



where E is the self-energy of <j). To the one-loop order, E 
is given by 

E w 2Alnf 6 . 



III. THE BOSE-BOSE MIXTURE 



The dilute BB mixture can be described by the action 



S 



where 



f P dTjd d a 



2m c 



(27) 



u =\Y1 ^i*«i 4 +«i2i^ii 2 i* 2 | 2 . 

a=l,2 

Here ^ a are the fields describing the bosonic atoms 
with masses m a . The couplings u a p are related 
to the scattering lengths a a p through the relations 
u±2 = 2uj±ai2 and u aa = 2ujj_a aa for the ID trap, 
and 1*12 = V / 47r(mi + m2)ui j_ / 'mim 2 ai 2 and u aa = 
yj 8ituj± /m a a aa for a 2D trap, where uj± is the trans- 
verse confining frequency. We shall take u a p > 0. 

By integrating out the fast modes, the one-loop RG 
equations in the dilute limit are completely identical to 
Eqs. (0| and ©, with the initial conditions 

r a (0) = (mi + m 2 )n a al , A Q(3 (0) = K d m af} u afi a . 

Here m a p = 2m a mp / {m\ + m 2 )- Hence, the fixed-point 
structure is identical to the one for the BF mixture. The 
solution of the RG equation for X a p{l) in d — 1 is of the 
form 



X a p(l) 



A* 



(28) 



6 



where 



A* 



Ca/3 



1 



Kim a paou a p 
while in d = 2 it is given by 

K 2 m a pu al 3 



1 + K 2 m a pu a pl 



(29) 



Scaling stops when 1 = 1*, where ri(/*) + r 2 (l*) = 1, 
yielding 



(30) 



(mi + m 2 )(^i + /U2)oo 



By a procedure similar to the Bose-Fermi mixture, £* in 
d = 1 can be expressed as 



mi 2 /(7rA») 



ao[(mi + 3m 2 )ni + (3mi + m 2 )n 2 ] 



(31) 



in the dilute limit. On the other hand, in d = 2, I* takes 
the form 



U = hi - 



toi 2 /(2-7t) 



[(mi + 3m 2 )ni + (3mi + m 2 )n 2 }a 2 



(32) 



in the extremely dilute limit. Here the extremely dilute 
limit is defined by I* ^ 1. (The diluteness condition only 
guarantees that e l * ^> 1.) 

We first use the above results to study the ID uniform 
mixture. Following the argument similar to the dilute 
BF mixture, the stability condition for the uniform mix- 
ture can be obtained from a mean-field analysis on the 
renormalizcd Hamiltonian, yielding 



4mim 2 
(mi + m 2 ) 2 



(33) 



where A*o = \ a p(l*)- It is clear that the fixed-point 
Hamiltonian (with A Q( g = A*) cannot satisfy Eq. (|33| . 
This means that the fixed-point Hamiltonian does not 
describe a uniform mixture, and thus we cannot take the 
continuum limit directly, as we have done for the BF 
mixture. 

To proceed, we insert Eq. (f28|) into Eq. (|33|) and use 
Eq. ([51]). Then, in the dilute limit (e u > 1), the stability 
condition can be expressed as 



< ma < U^ 



(34) 



where 
U^ 



Y 1 mmia (mi-roa) 2 iie 

27r 2 (mi +ni2)' 2 [(mi+3m2)ni + (3mi+m2)n2] 



is the critical strength of «i 2 in d = 1, and 
2mi 2 unu 22 



In the above, we have set A* = e = 1. From the expres- 
sion for Z4 (1) , we see that for the mixtures consisting of 
two different kinds of bosons, i.e. mi 7^ m 2 , the critical 
strength Uc is a linear function of boson densities 



U^ 



m 2 J 



mim 2 (mi — m 2 ) : 



■[(mi + 3m 2 )ni + (3mi + m 2 )n 2 ] 



in the dilute limit, which may be much smaller than u c . 
Moreover, Uc is is insensitive to the values of u aa . On 
the other hand, for the mixtures composed of two differ- 
ent hyperfine states of the same boson, i.e. mi = m 2 , 
u9^ reduces to 



UP 



"11 + U22 



in the dilute limit, which is independent of the boson 
densities. 

The stability conditions we obtained are valid only in 
the strong-coupling regime. We notice that they are dif- 
ferent from the result obtained from bosonization, which 

gave Uc 1 ^ = 7r 2 -\/nini/(mim 2 ). The study of the ID 
BB mixture with the exchange symmetry, i.e. mi = m 2 
and M11 = u 22 ~ U , by the finite-size density-matrix RG 
(FSDMRG) indicates that the phase separation occurs 
for all boson densities when u\2 > UJ^ Our conclusion 
is consistent with this result when the system possesses 
the exchange symmetry. 

Next, we turn into the 2D case. Due to the logarithmic 
RG flow of the coupling constants, we have two situa- 
tions: (i) When e l * ^> 1 but Z* = 0(1), we may take A*^ 
by their initial values, i.e. A*^ « K2vn a pu a ^. Inserting 
these into Eq. (|3"3"|) gives rise to the mean-field result 



< U 12 < y/UnU 2 2 



(35) 



(ii) In the extremely dilute limit, i.e. l„ 1, we insert 
Eq. (J29J) into Eq. ((33j) . Then, the stability condition in 
this situation becomes 



< U12 < Z4 (2) 



(36) 



where 



U. 



(2) 



1 



(mi — m,2 ) 2 u c 
167r(mi+m2) 



miun + m 2 u 22 



and h is given by Eq. (|32| . We see that for the mix- 
tures composed of two different hyperfine states of the 
same boson, i.e. mi = m 2 , the critical strength is given 

(2) 

by U c = u ci which is independent of the boson densi- 
ties, similar to the ID case. On the other hands, for the 
mixtures consisting of two different kinds of bosons, i.e. 
mi 7^ m 2 , Uc 2 ^ ~ 167r(mi +m 2 )/[(mi — m 2 ) 2 l*]. That is, 

it is a function of boson densities, which is insensitive to 

(2) 

the values of u aa - However, the dependence of U c on 
the boson densities is much weaker than the ID case due 
to the presence of the logarithm. This is a characteristic 
of the 2D Bosc gas. 
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A schematic phase diagram for the ID and 2D BB 
mixtures with repulsions is shown in Fig. [H where we 
have taken ni = n/2 = ri2- For mi ^ mi, Lie is a 

(2) 

linear function of n at low density and Uc is a linear 
function of 1/ In [1/(777100)] at extremely low density with 

7 ? = 2 7r(l+m 2 "/mi) :i - 0n th<3 ° ther ha " n ^ Wc = " c = 

in the limit na^ — > for mi = m2. In both cases, U^ 1 

approaches the mean-field value u m f = \Ju\\U2i at high 

density. We notice that the action S [Eq. ([27]) ] possesses 

an SU(2) symmetry when mi = m2, /-ti = p.2, and un = 

U22 = Ui2- In this SU(2) symmetrical limit, 

and the phase boundary may coincide with the mean-field 

prediction. 



IV. CONCLUSIONS AND DISCUSSIONS 

We employ the RG to study the dilute BF and BB 
mixtures with repulsive interactions in low dimensions. 
For the ID BF mixture, we show that the uniform mix- 
ture is stable against de-mixing for any finite strength 
of repulsions. This conclusion is consistent with the pre- 
diction obtained from the Bcthc ansatz when the system 
can be solved exactly. The bosonization approach gave 
a stability condition of the mean-field type at large bo- 
son densities (the weak-coupling regime)^, while at low 
boson densities (the strong-coupling regime) the stability 
condition becomes < U12 < ^ 2 V 'n b n ~t / '(mjm 7) . This 
result is inconsistent with the exact solution, while in the 
solvable limit ours is consistent with the Bcthc-ansatz re- 
sult. We believe that this inconsistency arises from the 
fact that the relations between the Luttinger liquid pa- 
rameters and the short-distance parameters like g bb and 
g b f are not reliable in the dilute limit where the system 
goes into the strong-coupling regime. 

For the 2D BF mixture with repulsions in the ex- 
tremely dilute limit, the uniform mixture is also stable 
irrespective of the interaction strength. The mean-field 
result is expected to be valid in the moderate and high 
density regime. We also calculate the critical tempera- 
ture T c for the SF phase. By comparing with T c for the 
single-component dilute Bose gas, we find that the inclu- 
sion of fcrmions increases T c . That is, the presence of 
fermions enhances the long range SF ordering. 

For the ID BB mixture, the low-energy physics is con- 
trolled by a fixed point which does not support the uni- 
form mixture. This observation immediately leads to two 
consequences: (i) First of all, the equilibrium properties 
of the uniform mixture will not exhibit universal behav- 
iors even in the strong-coupling regime because we can- 
not take the continuum limit, (ii) Next, the region of 
the stable uniform mixture with mi ^ m-2 can be very 
narrow in the dilute limit. For the 2D BB mixture in 
the extremely dilute limit, the region of the stable uni- 
form mixture with mi ^ m 2 is also narrow because the 
interactions are strongly renormalized in this limit. The 
mean-field result can be applied only in the moderate and 



high density regime. 

A recent paper also studied the stability conditions 
for the low-dimensional dilute BB mixture in terms of 
the RG approach^. The main distinction between our 
results and the ones in Rcf. 16 lies at the ID case with 
mi 7^ m-2, where our predicted value of Lie can be much 
smaller than u c in the dilute limit (our u c is basically 
u( ld ^ in Ref. 16), while in one dimension with mi = m2 
and in two dimensions, both are qualitatively consistent 
with each other. The discrepancy may be originated from 
the different criteria of the stability: In Ref. 16, the 
criterion for the stability is the positive definiteness of 
the potential term in the action at any scale I < I*. In 
our opinion, the use of the positive definiteness of the 
potential term in the action amounts to some kind of 
mean-field treatment, and such a procedure is meaningful 
only when the quantum fluctuations have been properly 
taken into account. Hence, this criterion can be used only 
at the scale I = 1* where the correlation length becomes 
of O(l) in the unit of the short-distance cutoff at I = I*. 
This is the criterion which we employed in this work. 
The issue about whether or not the critical strength for 
the ID uniform BB mixture can be very small in the 
dilute limit should be answered by further numerical or 
experimental studies. 
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Appendix A: Mean-field theory of the uniform 
Bose-Fermi mixture in d dimensions 

By taking the ansatz p = \<f>\ 2 , The mean-field free 
energy density / for the action S [Eq. (JXJ) ] is of the form 

f = + — p d{d + 2) 9(M/ - 9b fP) 
x(vf-g bf p) d/2+1 , 

where Q(x) = 0, 1 for x < and x > 0, respectively. The 
value of p is determined by the minimum of /, which 
leads to the mean-field equation 

d f 6/ Q(Mj - 9bfP)[2m f (iJ,f - g b fp)} d/2 = p-b - gbbP ■ 

The uniform mixture corresponds to the solution of the 
mean- field equation with < p < pf/g b f. Hence, for the 
uniform mixture, the free energy density is of the form 

g bb 2 2K d (2m f ) d / 2 d/2+1 

f M = -PbP + —P d{d + 2) ^ f ~ 9bfP > ' 

(Al) 
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where pb is the solution of the equation 



energy density 



Kd9bf 



[2m f (p f ~ 9bfP)] d/2 = Pb - 9bbP 



satisfying the constraint < p < pf/gbf. The density of 
fermions pf can be obtained from /m through the ther- 
modynamical relation pf = —dfni/dpf, yielding 



p f = ^[2 mf (p f -g bf p)}^ 



(A3) 



Using the thermodynamical relation p b = —dfuj/dpb 
where pb denotes the boson density, one may find that 
p = Pb- In terms of this relation and Eqs. (|A2|) . (|A3|) 
to eliminate p a , one may express /m as a function of 
p a - Using Euler's relation = —P and the thermody- 
namical relation P = —dE/dV, one may obtain the total 



(A2) 

V 



9bb 2 

-^Pb 



gbfPbPf 



where 



9ff 



1 / d 



2m f \Kd 



l+2/d 

-^9ffP f 



2/d 



(A4) 



By considering small density fluctuations, we obtain from 
Eq. (|A4p the linear stability condition 



2 „ tyff 2/d- 

9bf < —T-9bbn/ 



(A5) 
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